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PREFACE 


Certain sampling problems arising in connection with the 
Commission's continuous waybill study call for the frequent use of 
@ continuous series of random digits. Although various tables of 
random numbers have been published from time to time, none of those 
available contained a sufficient number of digits to serve the sam 
pling requirements of the waybill work. A newly developed method, 
outlined herein, was used to meet this need. 

The fact that random selection is a fundemental require- 
ment of any scientific sampling plan is generally recognized. Some 
arian failures in recent forecasts, due in part to imperfect sam- 

¥ pling procedure, have reemphasized the importance of sound sampling 
methods in statistical analysis. Tnasmach as the random numbers pro- 
duced for the purposes of the Commission's waybdill studies should be 
useful to statisticians generally, 105,000 of these numbers are pre- 
sented in this publication, together with a vrief technical discus- 
sion. 

The methodology used was devised by H. Burke Horten and 
R. fynes Smith III, and the tabulations were produced by the waybill 


unit of this Bureau under their supervision. 


Ww. H 8. Stevens, 


Director. 


TABLE OF 105,000 RANDOM DECIMAL DIGITS® 


Introduction 


The use of tables of random numbers in sampling and in other 
statistical problems has been frequently discussed, 1/ and tables of 
random numbers such as those compiled by Tippett, and by Kendall and 
Smith have been in use for some time. 2/ ‘The table of 105,000 numbers 
which follows may be used in the manner described in these various 
statistical discussions, the digits included having been developed by 
@ process which closely approximates an underlying rectangular distri- 
bution. It is also possible by the use of suitable transformations to 
obtain from the 105,000 numbers in this set, other sets of random num 
bers, which will follow any given distribution function. a/ 


The Randomizing Process 


: The mathematical principle involved in producing the accom 
panying set of random numbers has been explained elsewhere. Brief- 
ly, this process consists of adding groups of digits modulo their base 
(n) to obtain random digits for the final table. b/ By this method, 
the probability of each digit converges to 1/n under very general con- 
ditions. Thus, it is unnecessary to attempt to design a machine or a 
process which will have equal probabilities for each digit in single 
trials. Rather it is sufficient that the trials be independent and 
that the device (or process) be irregular in the sense that the proba- 
bility of obtaining each digit in each trial should be different fron 
zero, or more precisely, 


Pre2€ for all r and all s, 


where Prs represents the probability of obtaining the rtD digit in the 
eth trial, and € is a fixed positive number, arbitrarily small. c/ If 





* WNumbdered references appear at the end of the text. 
ee 
a/ Let y=/ ¢(z)dz. 
-© 


If y is distributed between 0 and 1 with uniform probability, x will 
be distributed according to #(x) , where # is any probability den- 
sity function. 


b/ This operation amounts to reading only the units position of a sum, 
For example, (34+74+9+4)mod 1023, this being the remainder after 
division of the sum by 10. 


c/ Actually, convergence to equal probabilities occurs under more gen- 
eral conditions, For details see first reference in 3/. 


this condition is satisfied, the probabilities of individual digits in 
the final set can be made as nearly equal as desired. his can be ac- 
complished by increasing the number of digits (original set) added 
modulo n to obtain each digit of the final set. 


Production of a Bandom Set 


In producing this particular set of random decimal digits it 
was calculated that addition modulo 10 of 10 decimal digits from the | 
original set would yield a set with probabilities in the range 1/10+10-”7 
Tis derived set was produced by tabulating 75,000 machine cards which 
had been punched from waybills received during three months in the 
regular course of the Commission's work. Fourteen colums, which ap- 
parently were independent, were selected from these cards. These in- 
cluded positions from such fields as the shipment weight, revenue, 
serial number of the car, etc. Bach of these columns was then tabu- 
lated by wiring into counters with no carryover and so controlled that 
the total in the units position of the counter would be cut in a sum 
mary card at the end of each 10 cards. ‘The totals (mod 10) of the 14 
columns for each block of ten cards were cut in the first 14 columns 
of the summary cards until 1,500 summary cards had been completed from 
the first 15,000 cards of the original set. The 1,500 summary cards 
were then replaced in the summary punch which was wired to punch in 
the next 14 columns, using totals from the second 15,000 cards of the 
original set. This operation was repeated until 70 columns of the sum 
mary cards had been punched. fhe tabulation following this discussion 
is a listing of the completed summary cards. 





fests for Randomess 


Any number of tests for local randomness may be devised. This 
particular set was subjected to three of the four tests used by Kendall 
and Smith. 1/ These tests, with the results of tests of significance 
are presented below. Table 1 shows the results of the frequency test in 
which the digits are merely counted and the totals ‘compared with the 
theoretical expectation of 10,500 of each. 


Table i 





Digit 0 1 2 3 4 5 


Actual frequency 10,420 10,496 10,540 10,507 10,445 10,569 
Expected f::squency 10,500 10,500 10,500 10,500 10,500 10,500 











6 7 8 9 Total 
Actual frequency 10,544 10,476 10,533 10,470 105,000 
Bxrpected freauency 10,500 10,500 10,500 10,500 105,000 





Test Results: xX*° <=1.938 P= .99 
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In the serial test adjacent pairs of digits are tabulated. 
According to the hypothesis of randomess, the frequencies of digits 
following any specified.digit should be equal, within tolerable sta- 
tistical limits. ‘the results of this test are presented in TPaple 2. 


Table 2 

4th. (t + 1) digit 
ligit 

0 1 2 4 5 6 7 g 9 Total 
0 O41, 1,073 1,005 1,095 1,016 1,021 1,017 1,063. 1,042 1,047 10,420 
1 1,030 1,024 1,043 1,024 1,057 1,081 1,117 1,043 1,037 1,040 10,496 
2 1,063 995 1,045 1,244 1,066 1,07% 1,005 1,063 1,056 1,032 10,5 
H 1,035 1,08 1,056 1,026 1,069 1,026 1,072 1,035 1,039 1,067 1 aa 

1 OR 1,061 1,106 1,028 1,011 1,082 1,028 983 1,053 1,058 10, 
5 1,049 1,065 1,08 1,046. 1,057 1,068 1,056 4058 1,060 1,0 205 oi 
6 1,046 1,060 1,070 1,102 1,038 1,029 1,063 1,080 1,008 1,048 
7 1,000 1,049 1,036 1,023 1,062 1,066. 1,080 1,053 1,079 1,028 0. 
& 1,091 1,078 1,004 1,027 1,019 1,032 1,065 1,071 1,088 1,058 10, ° 
9 1,030 1,009 1,089 992 1,050 1,093 1,042 1,047 1,071 1,048 10,470 


‘otal 10,420 10,496 10,540 10,507 10,445 10,569 10,544 10,476 10,533 10,470 105,000 





Test Results: x*- 75-461; P= .96 





In the gap test the number of non-zero digits appearing be- 
tween successive seros is tabulated. a/ The theoretical probability 
of occurrence of a gap of length k is given by p(k) = (9/10)1/10). 
Also, the probability of occurrence of a gap of length greater than 
or equal to k is given by: P(k) = (9/10)K It should be noted that 
two adjacent zeros are tabulated as a gap of length zero. ‘The total 
mumer of gaps between zeros is one less than the number of Zeros. 
The product of this total and the theoretical probabilities yields 
the expected frequencies for various gap lengths. The actual fre- 
quencies were compared with the theoretical frequencies, the results 
appearing in Table 3. The gep test could, of course, be applied to 
any digit, but due to the time required to make the test. only gaps 
between zeros were tabulated. 





a/ The gap test is one of a large group of tests which may be based on 
runs. 


Table 3 
Length of Gap O—] 2=3 4G 7-10 11-16 tT : Total 





Actual Frequency 1,962 1,620 1,857 1,685 1,505 1,792 10,419 
Expected Frequency 1,980 1,603 1,852 1,714 1,532 1,738 10,419 





fest Results: X* = 3.072; P= .69 





The machine methods used in running teots of local random 
ness did not lend themselves easily to tests of subsets of the 105,000 
digits. Tabulations for the serial and gap tests were run by colums, 
although the reading was horizontal. The tests of significance ap- 
plied to individual tabulations of the serial tests (by single columns) | 
and the gap tests (by groups of ten columns) gave a set of 77 proba- 
bilities which should be evenly distributed from 0 tol. A X2 test of 
the distribution of these probabilities is presented in Table 4. 


able 4 
Probability (tenths) 
Interval 
O-1 1-2 2- U5 5-6 6-7 7-8 8&9 9-10 Total 
Actual Frequency 8 9 5 9g § 10 7 24 %4 6 ai 


Expected Frequency 7.7 7-7 7.7 7-7 7-7 7-7 7-7 Teo? ToT 1? TT 


Test Results: y2 = 5.7273; P=.T7 


On the basis of the results of these tests it appears that the set of 
105,000 digits presented herein is locally randon. 


Use of Machine Methods of Selection 


In many sampling problems it is desirable to have sets of ran 
dom numbers in a form suitable for machine use. In most such cases ma- 
chines and sources of numerical data are available from which sets of 
any size can be constructed. Random selection of machine cards may alsc 
be accomplished without the use of recorded random numbers. Data re- 
corded on the cards themselves can be used to generate random digits by 
addition modulo 10, and these digits may be used to make the selection 


- FY = 


as the cards pass through the machine. For example, a listing of a 10 
percent random sample from a deck of cards could be obtained without 
disturbing the order of the cards. Several columns in the cards would 
be selected for independence and irregularity and the digits in these 
colums added into separate counters. The counters would not be clear- 
ed, so that at any time the digit in the units position of each counter 
would be the result of addition modulo 10 of the digits in the column 
wired to it. The units digits read successively from the counters 
would form a random sequence which could be used to control the listing. 
Extensions of this principle to other problems and for other types of 
equipment will snggest themselves. 


1/ Kendall, M, G,, and Smith, B. Babington, Randomness and random sam- 
pling numbers, Jour. Roy., Stat. Soc., Yol. 101 (1938), pp. 147 - 
166; Kendall,  G,, and Smith, B. Babington, Second paper on 
random sampling numbers, Jour. Roy., Stat. Soc. Supp., Vol. 6 
(1939), pp. 51 = 61. 


2/ Mippett, L. H. C., Random sampling numbers, Tracts for Computers, 

~ Wo. 15, (1927), Cambridges Kendall, M, G., and Smith, B. Babing- 
ton, Tables of random sampling numbers, Tracts for Computers, No. 
24 (1939), Cambridge. 


3/ H. Burke Horton and R. Tynes Smith III, A direct method for produc- 
ing random digits in any number system, Annals of Math. Stat., Vol. 
20 (1949), pp. 82 = 90; H. Burke Horton, A method for obtaining 
random numbers, Annals of Math. Stat., Vol. 19 (1948), pp. 81 - 85. 
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